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We study the effect of resonance decays on intensity in- 
terferometry for heavy ion collisions. Collective expansion of 
the source leads to a dependence of the two-particle corre- 
lation function on the pair momentum K. This opens the 
possibility to reconstruct the dynamics of the source from 
the K-dependence of the measured "HBT radii". Here we 
address the question to what extent resonance decays can 
fake such a flow signal. Within a simple parametrization for 
the emission function we present a comprehensive analysis of 
the interplay of flow and resonance decays on the one- and 
two-particle spectra. We discuss in detail the non-Gaussian 
features of the correlation function introduced by long-lived 
resonances and the resulting problems in extracting meaning- 
ful HBT radii. We propose to define them in terms of the 
second order g-moments of the correlator C(q, K). We show 
that this yields a more reliable characterisation of the corre- 
lator in terms of its width and the correlation strength A than 
other commonly used fit procedures. The normalized fourth- 
order (/-moments (kurtosis) provide a quantitative measure for 
the non-Gaussian features of the correlator. At least for the 
class of models studied here, the kurtosis helps separating ef- 
fects from expansion flow and resonance decays, and provides 
the cleanest signal to distinguish between scenarios with and 
without transverse flow. 



I. INTRODUCTION 

The only known way to obtain direct experimental 
information on the space-time structure of the particle 
emitting source created in a relativistic nuclear collision is 
through two-particle intensity interferometry QJ^]. This 
information is therefore indispensable for an assessment 
of theoretical models which try to reconstruct the final 
state of the collision from the measured single particle 
spectra and particle multiplicity densities in momentum 
space. Reliable estimates of the source geometry at par- 
ticle freeze-out are crucial for an experimental proof that 
high energy heavy ion collisions can successfully gener- 
ate large volumes of matter with high energy density. 
Direct information from two-particle correlations on the 
expansion dynamics at freeze-out further provides essen- 
tial constraints for theoretical models which extrapolate 
back in time towards the initial stages of the collision in 
order to make statements about a possible transition to 
deconfined quark matter. 

An important insight from recent theoretical research 



on Hanbury-Brown-Twiss (HBT) interferometry is that 
for dynamical sources which undergo collective expan- 
sion the HBT radius parameters, which characterize the 
width of the two-particle correlation function, develop a 
dependence on the pair momentum p|-|l3[]. The detailed 
momentum dependence is somewhat model-dependent, 
and in general it is not simple pT| . Still, it opens 
the crucial possibility to extract dynamical information 
on the source from interferometry data. Unfortunately, 
the most abundant candidates for interferometry studies, 
charged pions, are strongly contaminated by decay prod- 
ucts from unstable resonances some of which only decay 
long after hadron freeze-out [fbtpfl- Such resonance de- 
cays were shown to introduce an additional momentum 
dependence of the HBT radius parameters and of the in- 
tercept parameter ||,[l5| which complicates the extraction 
of the expansion flow. 

A systematic approach towards extracting the expan- 
sion velocity from experimental HBT data thus presup- 
poses a careful analysis of the interplay of flow and res- 
onance decays on the gross features of the two-particle 
correlation function. This is the aim of the present paper. 
We will use for our analysis a simple analytical model for 
the source function, which assumes local thermalization 
at freeze-out and produces hadronic resonances by ther- 
mal excitation. The model incorporates longitudinal and 
transverse expansion as well as a finite duration of parti- 
cle emission. The two most important parameters for our 
considerations, the temperature and transverse expan- 
sion velocity at freeze-out, can be varied independently. 
Our study thus complements published HBT analyses of 
source functions generated by hydrodynamic simulations 
where freeze-out is implemented along a sharp hypersur- 
face H and which do not allow easily to gain physical 
intuition by a systematic variation of the model param- 
eters. After freeze-out the resonances are allowed to de- 
cay according to an exponential proper time distribution 
along their trajectories, and the resulting emission func- 
tions of daughter particles (pions, kaons, etc.) are added 
to the direct emission function of particles of the same 
kind before calculating the correlation function. A dis- 
cussion of the momentum dependence of resonance de- 
cay effects on the 1- and 2-particle spectra requires the 
correct treatment of the decay phase-space [p[ p^Jl7| and 
does not permit the simplifying approximations leading, 
e.g., to Eq. (1) in Ref. @. 

The paper is organized as follows: In Sec. |TJ we review 
the extraction of space-time information on the source 
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from Gaussian fits to the correlation function. This cal- 



culational scheme is then extended in Sec. Ill to include 



sion function is sufficiently Gaussian one can approx- 
imate 



resonance decay contributions. The next three sections 
are devoted to a detailed model study based on this for- 
malism. 



C(q, K) « 1 + exp [-q^q„ (x»x v ) (K)} 



(2.4) 



In Sec. |V| we describe the model for the emis- where x^iK) = — 



sion function including resonance contributions. Results 
for the one- and two-particle spectra are presented in 
Sec. |y|, and a general discussion of the effects from res- 
onance decays on the shape of the correlation function 
is given there. In Sec. VI we then discuss in detail the 



x^{K) is the distance 



practical difficulties posed by the non-Gaussian features 
in the correlation function due to long-lived resonances, 
by co mpar ing different fitting procedures. This leads us 
in Sec. VII to the alternative method of ^-moments which 
provide a clean definition of the HBT radii and intercept 
parameter even for non-Gaussian correlation functions. 
These HBT radii show much weaker resonance decay ef- 
fects than the ones obtained in Ref. || by fitting a Gaus- 
sian function to a non-Gaussian correlator. The normal- 
ized fourth g-moment of the correlator, the kurtosis, pro- 
vides a quantitative measure for the deviations from a 
Gaussian shape as e.g. induced by resonance decays. We 
will show that, at least within the general class of source 
models studied here, the simultaneous study of the pair 
momentum dependence of the HBT radii, the intercept 
parameter and the kurtosis allows for a relatively clean 
separation of flow and resonance decay effects. We sum- 
marize our findings in Sec. VIII . The Appendix contains 



some background for readers interested in the technical 
details. 



II. GAUSSIAN PARAMETRIZATIONS OF THE 
CORRELATION FUNCTION 

For a given model for the emission function S(x,p) 
and assuming incoherent particle production as well as 
plane wave propagation, the invariant momentum spec- 
trum and two-particle HBT correlation function are given 
by f|§ 



Ep ~¥p = I (!> ' r Sir -i' 1 

C(q,K)f 



|J d A xS{x,K)e l " x \ 2 
\Jd 4 xS(x,K)\ 2 

l + \(e iq - x )\ 2 , 

/ d 4 x.f(x)S(x,K) 



(f(x))={f(x))(K) 



Jd 4 xS(x,K) 



(2.1) 

(2.2) 
(2.3) 



Eq. (^]^) is written down for identical bosons, and q = 
pi — p 2 , K = |(pi + p 2 ), with pi, p 2 on-shell such 
that K-q = 0. (f(x)) = (f(x))(K) denotes the (in- 
dependent) average of an arbitrary space-time function 
with the emission function S(x, K). As long as the emis- 



to the point x(K) of maximum emissivity of particles 
with momentum K in the source (the so-called "saddle 
point" of the source for particles with momentum K). 
In this approximation the two-particle correlation func- 
tion is completely determined by its Gaussian widths 
which in turn are directly given by the (if-dependent) 
space-time variances (x^x") of the emission function. 
The latter define the size of regions of homogeneity in 
the source [5| JlO|pl| , p^ | which effectively contribute to the 
Bose-Einstein correlations. Finer space-time structures 
of the source show up in deviations of the correlator from 
a Gaussian shape. 

In previous studies of analytically given emission func- 
tions, the correlator was sufficiently Gaussian to base all 
investigations on (2.4). Then one proceeds as follows: 
Due to the on-shell constraint K-q = only three of the 
four components of q which appear in the exponent are 
independent. The dependent component must be elimi- 
nated using the relation 



q° = (3 ■ q = 0±q o + ftqi . 



(2.5) 



Here /3 = K/K° » K/E K , with E K = y/m 2 +K 2 , 
is approximately the velocity of the pair, and we used 
the convention that I denotes the "longitudinal" (beam) 
direction (z-axis), o denotes the orthogonal "outward" 
direction (x-axis) which is oriented such that K = 
(K±, , 0, Ki) lies in the x-z-plane. Correspondingly (3 
has no y-component in the third Cartesian direction, the 
"sideward" dir ection: (3 S = 0. Due to the mass-shell 
constraint fl2.5]) , the inverse of the Fourier transform in 
( |2.2|) is not unique. The missing information required for 
the reconstruction of the (Gaussian) source in space-time 
from the measurable (Gaussian) HBT radii must thus be 
provided by model assumptions. 

In this paper we will deal only with azimuthally sym- 
metric sources for which the correlation function is sym- 
metric under q s — > — q s jp3|. Specifically, we will discuss 
two Gaussian parametrizations of C: 

1. The Cartesian parametrization H] is obtained by 



using (2.5) to eliminate q° in (2.4): 
C(q, K) = 1 + A exp [~q 2 s R 2 S (K) - q 2 Q R 2 (K) 
-2q oqi R 2 ol (K)] . 



(2.6) 



The corresponding size parameters are given by the 
space-time variances p4,pl| : 



R 2 S (K) 
R 2 (K) 

Rim 
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(2.7a) 
(2.7b) 
(2.7c) 
(2.7d) 
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For a detailed discussion of the meaning of these standard 
HBT parameters, in particular of the out-longitudinal 
{ol) cross term pjL and how they mix spatial and tem- 
poral aspects of the source, see Refs. ]io| , pd| , p3| . 

2. If one eliminates in (2.4) q Q and q s in terms of q±_ = 



can lead to appreciable non-Gaussian features in the cor- 



y/ q„ + Qg, q°, and qi one arrives at the Yano-Koonin- 



Podgoretskii (YKP) parametrization 



C(q, K) = 1 + A exp -i?i(K) q 



Rl(K){q?- q ° 2 ) 



(^o(K) 



i?|(K))(g.f/(K)r 



(2.8) 



where U(K) is a (i-T-dependent) 4- velocity with only a 
longitudinal spatial component: 



U(K) =7(K) (1,0,0, «(K)), with 7 
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(2.9) 



The YKP parameters R 2 ± (K), R$(K), and i?|(K) ex- 
tracted from such a fit do not depend on the longitudinal 
velocity of the observer system in which the correlation 
function is measured. They can again be expressed in 
terms of the space-time variances (x^x v ) and take 
their simplest form in the frame where w(K) vanishes 
p3| , ^2|p3| (the approximation in the last two expressions 
are discussed in (2||k|): 



Ri(K)=Ri(K) = (tf) 



Rj{K) 



^z-(0t/p x )x) 
Rl{K) = ((t-x/^) 2 )-(f)/Pl 



-{Pi/0±) 2 (f) 
(t 2 ). 



(2.10a) 

(2.10b) 
(2.10c) 



The expressions (2.7,2.10) for the HBT parameters are 
useful for two reasons: (i) They result in an apprecia- 
ble technic al si mplification because instead of the Fourier 
transform ( 2f2 ) only a small number of 4-dimensional real 
integrals over the source function must be evaluated to 
completely determine the correlation function. Their ac- 
curacy has been checked in (ll]] for models of the type 
to be used below and, in the absence of resonance de- 
cays, for hydrodynamic sources with a sharp freeze-out 
hypersurface in P8| . (ii) They provide an intuitive under- 
standing of which space-time features of the source are 
reflected by the various g-dependencies of the correlator. 
However, their range of validity is limited by the fact that 
strictly speaking the space-time variances determine only 
the curvature of the correlator at q = 0: 



(fa - 0ii)(xj - Pjt)) = - 



l5 2 C(q,K) 



2 dqi dqj 



(2.11) 



q=0 



This agrees with the widths of the correlator if and only 
if C(q, K) is Gaussian. We will see that resonance decays 



relation function. If this is the case, Eqs. (2.7) and (2.1C) 
no longer provide quantitatively reliable expressions for 
the observed half widths of the correlator. A detailed 
discussion will follow in Sec. [y|. 

III. RESONANCE DECAY CONTRIBUTIONS 



We concentrate on charged pion {tt + tt + or 7r~7r~) cor- 
relations. In the presence of resonance decays, the emis- 
sion function is the sum of a direct term plus one addi- 
tional term for each resonance decay channel with a pion 
of the desired charge in the final state: 

S„{x,p) = Sf T (x,p) + S r ^(x,p) . (3.1) 

r^TT 

Note that the sum is over decay channels, not just 
over resonances. We compute the emission functions 
S r ^ w (x,p) for the decay pions from the direct emission 
functions S^ n '{X,P) for the resonances taking into ac- 
count the correct decay kinematics for 2- and 3-body de- 
cays, 

5 r _»(x,p) =J2 I I drTe- TT S^ (x~^tJ" 



R JO 



M 



(3.2) 



From now on capital letters denote variables associated 
with the parent resonance, while lowercase letters denote 
pion variables. Here, T is the total decay width of the res- 
onance, and J R goes over the kinematically allowed 
resonance momenta as described in Appendix |A|. Please 
note that the momenta p and in this expression are in 
general different, in contrast to the approximation used 
in Eq. (1) of Ref. [p^| . This is important for the following 
discussion of the momentum dependence of the correla- 
tor. 

The complete two-particle correlation function is then 
given by 



C(q,K)=l 



+2 



\Sf r {q,K)f 
1^(0, K)\* 

J2r^Ms^{q,K)s r ^{q,K)} 



+ 



\S„{0,K)\2 
M,K)\ 2 



(3.3) 



1^(0, K )|2 

where the denominator includes all resonance contribu- 



tions according to (3.1). The last term in the numerator 
can be neglected if resonance production is small [ B9[ . 
However, in ultrarelativistic heavy ion collisions a major 
fraction of all final state pions stem from resonance de- 
cays (see Fig. [I]) and this "Grassberger approximation" 
cannot be used. 
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For later reference, we extend the expressions given 
in Sec. || for the HBT parameters in terms of space- 
time variances of the source to include resonance decay 
contributions: 

Here the sum runs over all contributions, including the 
direct pions. It is instructive to rewrite the average over 
the emission function in the following form: 



{x v ){K)=Y,U{K) WrW, 

r 

(x^)(K) = fr(K) { Xlx X v ) r (K) , 
r 

where we introduced the single-particle fractions j 
/ d 4 xS r ^{x,K) _ dN^/d 3 K 



(3.5) 



fr(K) = 



J2 r J d±x S r ^(x, K) dN™/d 3 K 

E. 



(3.6) 



These give the fraction of single pions with momentum 
K resulting from decay channel r. We also defined the 
average (. . .) r with the effective pion emission function 
arising from this particular channel: 



{...)r{K) = 

The variances (m 



J d 4 x . . . S r -nr(x, K) 
J d 4 xS r ^ 7r {x,K) ' 

can then be rewritten as 



(3.7) 



(x^Xy) — ^ f r (x ^Xij) r -f~ ^ fr {3'r.r' fr' ) { 



(3.8) 



The first term has a simple intuitive interpretation: each 
resonance decay channel r contributes an effective emis- 
sion function SV->7r- The full variance is calculated by 
weighting the variance (homogeneity length) of the emis- 
sion function from a particular decay channel r with the 
fraction f r with which this channel contributes to the 
single particle spectrum. However, the different effective 
emission functions S r -^(x,p) have in general different 
saddle points; this gives rise to the second term in (3.5) 
which somewhat spoils its intuitive interpretation. 

Also, the full emission function (3.1) is a superposition 
of sources with widely differing sizes since long-lived res- 
onances contribute long exponential tails to the emission 
function Sy—^r p|,p^|. It is easy to see that this leads to 
non-Gaussian correlation functions: Consider a simple 1- 
dimensional toy model where the emission function is a 
sum of two Gaussian terms, one of width i?dir for direct 
pions and one of width i?haio for pions from a resonance, 
with weights 1 — e and e, respectively: 



S n (x, K) = S%*{x, K) + S r ^(x, K) 

= (1 - e) e -* 2 /(2flD + e e -* 2 /(2fliU) . (3.9) 



According to ( |3.3| ) the correlator is then a superposition 
of three Gaussians which for i?haio S> -Rdh- have very 
different widths: 



2 e -K„<t 



C(q,K)-l = (l-e 

+2e(l-e)e 



e 2 e -K &l y 



-(«L+-RLlo)9 2 /2 



(3.10) 



Obviously, if e is small, the rough structure of the cor- 
relator will be determined by the large and broad direct 
contribution. The two other contributions will, however, 
modify its functional form as follows: 

(i) If the resonance is shortlived such that i?haio ~ -Rdir, 
its effect on the correlator will be minor; its shape will 
remain roughly Gaussian, with a width somewhere be- 
tween l/i?dir and l/i?haio, depending on the weight e of 
the resonance contribution. 

(ii) If the resonance lifetime and thus i?i ia io are extremely 
large, the second and third term in ( 3. 1C ) will be very nar- 
row and, due to the finite two-track resolution of every 
experiment, may escape detection; then the correlator 
looks again Gaussian with a width l/i?dir, but at q = 
it will not approach the value 2, but 1 + (1 — e) 2 < 2. The 
correlation appears to be incomplete, with a "correlation 
strength" A = (1 — f r ) 2 = (1 - e) 2 . 

(iii) If the resonance lifetime is in between such that 
-Rhaio ^S> i?dir but 1 /i?haio is still large enough to be 
experimentally resolved, all three Gaussians contribute, 
and the full correlator deviates strongly from a single 
Gaussian. 




500 1000 1500 
K ± (MeV) 



1.0 2.0 3.0 

y 



thermal p A K L m t\ r\' K° L, A 
FI G. 1 . The resonance fractions fr(y,p±) according to 
Eq. (3.6) for T = 150 MeV. Upper row: no transverse flow, 
r\j = 0; lower row: r\j = 0.3. Left column: f r as a function 
of transverse momentum at central rapidity; right column: f r 
as function of rapidity at p± = 0. 
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In cases (ii) and (iii) the space-time variances calcu- 
lated from Eq. yield misleading or outright wrong 
results for the width of the correlation function. They 
give the curvature of the correlator at q = 0: 



/dir(£ 2 )dir + /halo (x 2 ) halo — (1 — e)-Rdir + e ^halo 



1 d 2 C(q) 
2 



dq 2 



(3.11) 



g=0 



In case (ii), for not too small values of e, this is dom- 
inated by the second term although the resonance con- 
tribution is not even visible in the measured correlator! 
On a quantitative level, the situation is not very much 



better for case (iii) (see Sec. VI A 1 for a more detailed 
discussion) . 

However, if the resonances can be clearly separated 
into two distinct classes, one with very short lifetimes 
of order 1 fm/c, the other with very long lifetimes of 
order 100 fm/c or longer, then space-time variances can 
again provide an accurate measure for the width of the 
correlation function. To achieve this, one must leave out 



the long-lived resonances from the sum over r in (3.8), i.e. 
one restricts the calculation of the space-time variances 
to the "core" of the emission function from direct pions 
and short-lived resonances J|(]|5l| . Since the contribution 
from long-lived resonances to the correlator cannot be 
resolved experimentally (while they do contribute to the 
single-particle spectra), one includes them via a reduced 
correlation strength A: 



X(K)=[l- Yl fr(K) 



(3.12) 



r=longlivcd 



The ^-dependence of A will be discussed in Sec. |y|. 

The real problem comes from resonances with an inter- 
mediate lifetime. They cause appreciable deviations from 
a Gaussian behaviour for the correlator and cannot be re- 
liably treated by the method of space-time variances. In 
nature there is only one such resonance, the u> meson 
with its 23.4 fm/c lifetime. At low K± it contributes up 
to 10% of all pions (/ W (K = 0) w 0.1), and their non- 
Gaussian effects on the correlator can be clearly seen. 



They will be discussed extensively in Sections VI , VII 



IV. 



A SIMPLE MODEL FOR THE EMISSION 
FUNCTION 



As discussed after Eq. (2.5), a completely model- 
independent HBT analysis is not possible. In this Section 
we define a simple model for the emission function in rel- 
ativistic nuclear collisions which will be used in the rest 
of the paper for quantitative studies. It has been used 
extensively in the literature [p|- |l3| , p3"| | , and we present 
a simple extension to include resonance production. It 



implements the essential features expected from sources 
created in nuclear collisions: It assumes local thermal- 
ization prior to freeze-out and incorporates its collective 
expansion in the longitudinal and transverse directions. 
On the geometric side, the source has a finite size in the 
spatial and temporal directions, i.e. it implements a fi- 
nite, but non-zero duration for particle emission. 

The emission function for particle species r is taken as 



S^(x,P) 



where 



H(x) 



7t(At) 



^MxcoshCT-,) 
x exp 

exp 



P ■ u(x) — IA r 
T 



H(x) (4.1) 



r 2 _ (r? - 7?o) 2 _ (t - Tpf 
2R 2 2(Ar?) 2 2(At) 2 



(4.2) 



with proper time r = \fW 
V = £log[(f- 



z 2 and space-time rapidity 
z)/(t— z)]. The physical meaning of the 
parameters has been explained in detail in Refs. [p|- p3p3 
to which we refer the reader. The only new ingredi- 
ents are a factor 2J r + 1 for the spin degeneracy (due 
to charge identification in the experiment each isospin 
state must be treated separately) , and a chemical poten- 
tial [i r for each resonance r. This means that all particles 
are assumed to freeze out with the same geometric char- 
acteristics and the same collective flow, superimposed by 
thermal motion with the same temperature. The pos- 
sible consequences of particle-specific freeze-out j$4 35 
will have to be discussed elsewhere. 

For later reference we note that the function H(x) is 
normalized to the total comoving 3-volume according to 



S d ' xH{x) 



Vrms i 



us rms ; 

r? rms = At] . 



(4.3) 
(4.4) 



Note that the rms widths in x- and y-direction are each 
given by R. If the Gaussians in H{x) were replaced 
by box functions ]3^ , |33[ |, the equivalent box dimensions 
(with the same rms radii) would be R = 2R, fj — V3 Ar]. 

For the flow profile we assume Jl2j Bjorken scaling in 
the longitudinal direction, vi = z/t, and a linear trans- 
verse flow rapidity profile p6|: 



Vt(r) = rif — 



(4.5) 



In spite of the longitudinal boost-invariance of the flow, 
the source as a whole is not boost-invariant due to the 
finite extension in r\ provided by the second Gaussian in 



(4.1) 



Inserting the parametrization (A5) for P the emission 
function (4.1) becomes pl|] 
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2 7-1-1 

Sf T (x, P) = r T" M T cosh(F - rf) e^ /T H(x) 



(2tt) 3 



/i r = b r /J, B + SrHS 



(4.7) 



x exp 



Mi 



cosh(y — rf) coshr)t(r) 



Pt 



+ — smhr) t (r) cos((f> - $)J . (4.6) 

The direct pion component S^"(x,p) is obtained from 
this expression by setting r = n, P = p, = 0, [i w = 0, 
and $ = (see ( A4)). This last condition reflects a 
choice for the orientation of the coordinate system such 
that the transverse momentum pj_ of the decay pion lies 
in the x — z plane. For the transverse momentum 
of resonances which contribute pions with the same 
as the directly emitted ones, in general a non-vanishing 
azimuthal angle is required, see Appendix [A|. 

The chemical potentials [i r will be fixed by the as- 
sumption of chemical equilibrium at freeze-out. Then 
baryon number and strangeness conservation in the fire- 
ball demand the existence of two independent chemical 
potentials [Ib and fis> with 



Decay channel r 


M (MeV) 


T (MeV) 


J 


6 r-„r- 


P 




TV 7T U 


770 


150 


1 


1.0 


P° 




7T~7T + 


770 


150 


1 


1.0 


A' 




> ir~n 


1232 


115 


3/2 


1.0 


A° - 


► 7T~p 


1232 


115 


3/2 


(1/3) x 1.0 


A 4 




> 7r~n 


1232 


115 


3/2 


(1/3) x 1.0 


A+ 


+- _ 




1232 


115 


3/2 


1.0 


K*° 


— > 




892 


50 


1 


(2/3) x 1.0 


K*~ 






892 


50 


1 


(2/3) x 1.0 


£*" - 


• 7T 


~A(1116) 


1385 


36 


1/2 


0.88 


E*~ -> 




£°(1193) 


1385 


36 


1/2 


(1/2) x 0.12 




7r~ 


E+(1193) 


1385 


36 


1/2 


(1/2) x 0.12 


E* + - 


' 7T 


-A(1116) 


1385 


36 


1/2 


0.88 


E* + - 




'S°(1193) 


1385 


36 


1/2 


(1/2) x 0.12 


E*° -» 


7T~ 


E"(1193) 


1385 


36 


1/2 


(1/2) x 0.12 


U) — 


IT 


-7v+n° 


782 


8.43 


1 


0.89 



i] — ► n n r\ 



E ->i n 
£+ -» 7r~n 
E° -> 7A -> p7T~ 

A —> pir~ 



547 

958 

498 

1193 
1193 
1193 
1116 



1.2 x 10" 

0.2 

« 

« 
« 
« 
« 







1/2 
1/2 
1/2 
1/2 



0.24 

0.44 

0.69 

1.0 
1.0 

0.65 
0.65 



TABLE I. The resonance decay contributions to n~ pro- 
duction considered in the present work. Where applicable the 
factor in front of the branching ratio is the Clebsch-Gordon 
coefficient for the particular decay channel. 



where b r and s r are the baryon number and strangeness 
of resonance r, respectively. The condition of overall 
strangeness neutrality of the fireball allows to eliminate 
[is in terms of T and [ j37| . 

Unless stated otherwise, the numerical calculations |44| ] 
below are done with the set of source parameters T = 150 
MeV, R = 5 fm, Arj = 1.2, r = 5 fm/c, At = 1 fm/c 
and /its = [is — 0. We will work in the fireball cm. 
system and thus set rjo = 0. 

The resonance channels included are listed in Table Q. 
The £(1193) and A(1116) are treated as one baryonic 
resonance F(1150) at an average mass of 1150 MeV. For 
simplicity the decay cascade E° — > 7A — > pir~ is re- 
placed by an effective two-particle decay S° — > j>7r~, since 
the photon in the S°-decay is known not to change the 
shape of the hyperon spectrum |j32|| . The n~ decay con- 
tributions from the cascades rj' — ► ... +77 — * 7r~ + ... and 
S* — ► ...+ Y(1150) — > 7r~ + ... are taken into account by 
enhanced branching ratios for the 1" and 77 decay chan- 
nels. These crude approximations are not problematic 
because they concern quantitatively small contributions. 
The cascade decays just mentioned affect the intercept 
parameter on the level of a few percent; the K- depen- 
dence of the HBT radius parameters remains essentially 
unaffected. K\ decays are neglected because the long 
K\ lifetime (cr = 15.5 m) makes them invisible for most 
detectors. 



V. RESULTS FOR ONE- AND TWO-PARTICLE 
SPECTRA 

We now present a quantitative analysis of the one- and 
two- particle spectra for the model described in Sec. IV. 
Both types of spectra can be expressed in terms of the 
four-dimensional Fourier transforms of the direct emis- 
sion functions, S^"(q, P^), see Appendix [a|. We show in 
Appendix [b| how the latter can be reduced analytically 
to 2-dimensional integrals over r and 77: 



i j^M^Jdr 1 (l + iA q ) 



x cosh(r; - Y) e lATa e 
r dr e 



At„-±A 2 (A 



V 



cosh(r/—Y) cosh r}t 



h (y/C-W± 



where 



C(r) 



P? 2 
= Tp-smh Vt{r) 

Pi 
T 



2 2 
r q x , 



(5.1) 



(5.2a) 



D ± (r) = -2r-^-smhr) t (r) 
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x (q Q cos$± + <7 s sin<f>±) 
A(rj) = (q cosh 77 — qi sinhry) , 

A q (r))=A(r,) < > AT) ~ 



TO 



(5.2b) 
(5.2c) 

(5. 2d) 



The Bessel function 1$ arises from the 0-integration 
while the terms containing A and A q stem from the r- 
integration. Please note that the azimuthal rotation of 
the resonance transverse momentum relative to the 
pion transverse momentum pj^ (which defines the x-axis 
of our coordinate system) enters only through the com- 
bination in brackets in Eq. ( 5.2b| ); the latter stems from 
the scalar product q ■ P^ , see Eq. ( |A19| ) . This means that 
the dependence on $± can be shifted from to qj^ by 
a common rotation by the angle $±: 

Sf r (q° , 9o , 9s , ft ; £ p , P ± cos $ ± , P ± sin *±,P L ) = 
Sf 1 ' (q ,q o cos $± + q s sin $±, q s cos $±- 
q o sm<S>±, qi ;E p ,P ± ,0,P L ) . (5.3) 

Note that this identity does not depend on the model 
for the emission function. It shows that the resonance 
decay kinematics leads to a mixing of the sideward and 
outward g-dependencies of the correlation functions that 
would be obtained from the resonances if one could use 
them for interferometry directly. This feature is lost in 
the approximation leading to Eq. (2) in Ref. p8| . 

For t he dir ect pion contribution, $± is to be set to zero 
in Q, (|J). 



A. The resonance fractions f r (K) 



The single particle momentum spectrum (2.1) is the 
space-time integral over the emission function ( |3.l|) , 



7r dy dm^ 



= / d xS n (x,p) = S n (q = 0;m±,y). (5.4) 



It is thus given by the Fourier transform (5.1) of the 
emission function at zero relative momentum. From this 
expression it is straightforw ard to evaluate the resonance 
fractions / r (j/,roi) of Eq. (3.6). For later referene they 
are shown in Fig. EL At central rapidity and small trans- 
verse momentum in our model only about 40% of the 
pions are emitted directly while more than half of the 
pions stem from resonance decays. The direct fraction 
increases rapidly with increasing transverse momentum, 
but very slowly with increasing longitudinal momentum 
resp. rapidity. In fact, most resonance fractions are 
nearly independent of rapidity || . At large p± the res- 
onance contributions to the single particle spectrum die 
out |L7]]. The largest resonance contribution comes from 
the p meson, due to its relatively small mass and large 
spin degeneray factor. The 7/, which is still lighter, has 
no spin and a small branching ratio into pions. As can 



be seen in the lower row of Fig. [j] the resonance fractions 
are only weakly affected by transverse flow: at small p± 
the direct fraction increases slightly while at large p± the 
tendency is opposite (see Sec. VB). 



B. Single particle transverse momentum spectra 



Integrating (5.4) over rapidity we obtain the single par- 
ticle transverse momentum distribution: 



dN^ 
dm\ 



tt J dyS« 1T (0;y,m ± ) 

+ ^7r J ' dyS r ^(0;y,mj_) . 



(5.5) 



The resonance decay contributions are given according 
to Eqs. jAl8|) and (|Al|) by 



S r ^(0: y, mj = 2M / S r dir (0; Y, M±) . 



(5.6) 



R 



The factor 2 results from the sum over <I>±, noting that 
at q = the integrand is independent of &± ( see Ap- 
pendix |b|). Writing Y = y + vAY (see Eq. ( A16 )) where 
Ay is independent of y, the y-integration can be pulled 
through the integrals J R over the decay phase space, 
yielding fT§ 



drn?^ 



dm 



dm i 



2M r 



dm 



R 



dMj 



(5.7) 



The transverse momentum spectra of the directly emitted 
resonances r are given by expression (B4) |l^,^2|: 



dm 



1 



dM\ 



At: 2 

x Mi 



(2ttR 2 ■ 2r Ar/) e 



Mr/T 



x/ (^sinhr ?t (C)) 



tfi^coshijttO) 
(5.8) 



where we substituted £ = r/R under the integral. Note 
that the geometric parameters R, Arj, tq of the source 
enter only in the normalization of the spectrum through 
the effective volume (4.3). Thus the shape of the (y- 
integrated!) single-particle transverse momentum spec- 
trum contains no information on the source geometry, in 
agreement with general arguments presented e.g. in p2[ . 
According to (pV? ,5.8), the unnormalized transverse mo- 
mentum dependence is fully determined by the rest mass 
M, the temperature T (or T(£) if T were r-dependent), 
and the transverse flow profile %(£,) = ?7/C™- 

For later reference we plot in Fig. |^ the pion trans- 
verse mass spectrum for the two sets of source parame- 
ters for which we compute two-particle correlations be- 
low. All resonance decay contributions are shown sep- 
arately. The only 3-body decays are those of the w, 77 
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and rj' whose decay pions are seen to be particularly con- 
centrated at small p±. (A similar low-p_L concentration 
occurs for pions from Kg decays, due to the small decay 
phase space in this particular 2-body decay.) Compar- 
ing the top panel (no transverse flow, rjf = 0) with the 
bottom panel (rjf = 0.3) one observes the well-known 
flattening of the transverse mass spectrum by transverse 
radial flow [{32|-|34 138 1 . The direct pions reflect essentially 



an effective "blueshifted" temperature T e g — r ^j^j 

[gSj. But the heavier resonances, in the region P± < M r , 
are affected much more strongly by transverse flow since 
at small P± the flattening of the spectra by flow is pro- 
portional to the particle rest mass p-JpSl . 




m ± - m K (MeV) 

FIG. 2. The single pion transverse mass spectrum for 
T — 150 MeV and jib = A*s — 0. The overall normaliza- 
tion is arbitrary, the relative normalizations of the various 
resonance contributions are fixed by the assumption of ther- 
mal and chemical equilibrium. Upper panel: no transverse 
flow, rjf = 0; lower panel: rjf = 0.3. 

Fig. ||b shows that this effect on the parent resonances 
is also reflected in the spectra of the daughter pions, ex- 
plaining the slight rise with rjf of the resonance fractions 
at large m ± . This flattening of the transverse mass spec- 
tra by transverse flow, suggested in Refs. [^8| , p2|j3^ ] as an 
explanation for the observed features of the single particle 
spectra from 28 Si- and 32 S-induced collisions at the AGS 
and SPS, seems to have been confirmed by recent colli- 
sion experiments with very heavy ions (Au+Au at the 



AGS, Pb+Pb at the SPS, see contributions by Y. Ak- 
iba, R. Lacasse, Nu Xu and P. Jones at the recent Quark 
Matter '96 conference (3^]). One of the main goals of 
two-particle interferometry is to obtain an independent 
and more direct measure of the transverse expansion ve- 
locity at freeze-out, to confirm this picture and further 
discriminate against possible alternative explanations. 

T|f = 0, Y = 

q s = q ,=o q =qi=o q =q s =° 



- 1.0 
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10 20 30 40 10 20 30 40 10 20 30 40 50 
q^MeV) q s (MeV) q|(MeV) 

FIG. 3. The two-particle correlator C(q, K) for ir~ pairs 
with pair rapidity Y = in the CMS. Each row of diagrams 
corresponds to a different value for the transverse pair mo- 
mentum K± (K± = 0, 200, 400, 600 and 800 MeV from top to 
bottom). Left column: the correlator in the outward direction 
at q s = qi — 0. Middle column: the correlator in the sideward 
direction at q — qi = 0. Right column: the correlator in the 
longitudinal direction at q 3 = q — 0. Source parameters as 
in Sec. |iy|, the transverse flow rjf has been set to zero. 
Here and in the following plots the different lines have 
the following meaning: Thin solid line: thermal pi- 
ons only. Long-dashed: including additionally p-decays. 
Short-dashed: including additionally all other shortlived res- 
onances (A,_K"*,E*, see Table Q). Dash-dotted: adding also 
w-decays. Thick solid line: adding also all longlived reso- 
nances (rj, rj' , Kg, E, A, see Table mj. 



C. Two-particle correlations 

In Figs. [| and |J we plot the two-pion correlator 
C(q, K) in the three Cartesian directions of q for zero 
and non-zero transverse flow rjf, respectively. We use 
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the letter Y to denote the rapidity of the pair, and K± 
for its transverse momentum (transverse mass). 
The pion pairs in Figs. || and [| have pair rapidity Y = 
in the CMS, and transverse momenta ranging from to 
800 MeV/c (top to bottom). The correlation functions 
were calculated by numerically evaluating Eq. (|3.3|) for 



the source parameters given in Sec. IV. 
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FIG. 4 
rif = 0.3. 
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Same as Fie 
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|j| but for nonzero transverse flow 



Within each pl ot, t he different lines show the effect 
of adding in Eq. (3.3) in the sum over decay channels 
r successively more resonances (see Table ||): first the 
abundant p, then the other short-lived resonances, then 
the u> with its intermediate lifetime, and finally all the 
long-lived resonances. Comparing these plots row by row 
gives one a feeling for the fG_-dependence of the correla- 
tion function and the various resonance contributions. In 
the following two subsections we give a rough and general 
discussion of the main features of the correlator without 
and with transverse flow of the source, respectively, be- 
fore proceeding to a quantitative analysis in Sec. VI. 



No transverse flow (Fig. 



The direct thermal contribution leads to a correlation 
function with a nearly Gaussian shape in all directions 
qi, i — o, s, I, and for all pair momenta K±. As Kj_ in- 



creases, the correlator becomes rapidly wider in the lon- 
gitudinal direction while in the two transverse directions 
the changes are hard to see and require a finer analy- 
sis (Sec. VI). As more and more of the short-lived res- 
onances are added, the width of the correlator becomes 
smaller, again with a larger effect in the longitudinal than 
in the two transverse directions. A much stronger effect 
is caused by the w-meson; now the narrowing of the cor- 
relation function is also clearly seen in the transverse 
directions, and the correlator becomes markedly non- 
Gaussian. As the long-lived resonances are added, the 
intercept A of the correlator at q = decreases below 1. 
This is a matter of q- resolution (we stop at |q| = 1 MeV) 
- the contribution from the long-lived resonances is en- 
tirely concentrated in a 5-function like structure near the 
origin, and with infinite resolution the correlator could 
be seen to actually reach the value 1 at q = 0. This is, of 
course, an extreme deviation from Gaussian behaviour. 

Long-lived resonances thus lead to apparently incom- 
plete correlations, A < I This effect becomes even 
stronger, if the correlator is projected onto one particu- 
lar g-direction by averaging over a finite window in the 
other directions where the correlator has already dropped 
below A @. 

As the pion pair momentum K± increases, all reso- 
nance effects on the width and strength of the correlator 
are seen to decrease. This is a direct consequence of the 
decreasing resonance fractions, see Fig. [l| 

The above lifetime hierarchy of resonance effects can 
be understood in terms of the following simple picture: 

• Short-lived resonances, T > 30 MeV. In the rest 
frame of the particle emitting fluid element these 
resonances decay very close to their production 
point, especially if they are heavy and have only 
small thermal velocities. This means that the emis- 
sion function S^- nr of the daughter pions has a very 
similar spatial structure as that of the parent reso- 
nance, S^ lr , although at a shifted momentum and 
shifted in time by the lifetime of the resonance. As 
only R and Ri are sensitive to the lifetime of the 
source, the shift in time affects the correlation func- 
tion only in the outward and longitudinal direc- 
tions. The stronger effect on Ri (which is obvious 
from the right column in Fig. ||) is a consequence 
of the boost-invariant longitudinal expansion of our 
source: as the decay pions are emitted at a later 
proper time t, and since the longitudinal length of 
homogeneity increases with r because the longitu- 
dinal velocity gradients decrease [|l0| , the decay pi- 
ons show a larger longitudinal homogeneity length 
than the direct pions. - Since the Fourier transform 
of the direct emission function is rather Gaussian 
and the decay pions from short-lived resonances ap- 
pear close to the emission point of the parent, they 
maintain the Gaussian features of the correlator. 
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• Long-lived resonances, r< 1 MeV. These are the 
77 and 77', with lifetimes err ~ 17.000 and 1000 
fm, respectively, and the weak decays of Kg and 
the hyperons which on average propagate several 
cm. (The decays of K\ and charged kaons are 
not included in our calculation because they are 
reconstructed in most experiments.) Even with 
thermal velocities these particles travel far outside 
the direct emission region before decaying, generat- 
ing a daughter pion emission function S r ^ n with a 
very large spatial support. The Fourier transform 
S r ^Tr(q, K) thus decays very rapidly for q ^ 0, giv- 
ing no contribution in the experimentally accessi- 
ble region q > 1 MeV. (This lower limit in q arises 
from the finite two-track resolution in the experi- 
ments.) The decay pions do, however, contribute 
to the single particle spectrum S r ^ 1T (q = 0,K) in 
the denominator and thus "dilute" the correlation. 
In this way long-lived resonances decrease the cor- 
relation strength A without, however, affecting the 
shape of the correlator where it can be measured. 

• Moderately long-lived resonances, 1 MeV < Y < 30 
MeV. There is only one such resonance, the u> me- 
son. It is not sufficiently long-lived to escape de- 
tection in the correlator, and thus it does not af- 
fect the intercept parameter A. Its lifetime is, how- 
ever, long enough to cause a long exponential tail 
in Sw-nr(x, K). This seriously distorts the shape of 
the correlator and destroys its Gaussian form. 



2. Non-zero transverse flow (Fig. 

The main difference between Figs. || and ^ is that the 
effects from the short-lived resonances and the u> on the 
shape of the correlator are weaker. The primary re ason 
for this behaviour is that for the class of models ( |4.6| ) 
the transverse size Rt of the effective emission region for 
heavy resonances shrinks for nonzero transverse flow. In 
Gaussian saddle point approximation, this transverse size 
Rt can be calculated from S^"(x,P) in (4.6) as 



Ri 



R 



+ (M ± /T) V j 



(5.9) 



This is not accurate enough for quantitative studies ]Tl| ] 
but gives the correct tendency and right order of mag- 
nitude. Going like 77^ this effect is small, but it tends 
to increase the width of the correlator, counteracting the 
basic tendency of resonance contributions to make the 
correlator narrower. For 77/ = 0.3 the two effects are 
seen to more or less balance each other in the sideward 
correlator, leaving practically no trace of the shortlived 
resonances including the uj. A similar effect is seen in 



the outward and longitudinal directions, but there the 
dominant lifetime effect discussed above prevails. 

Please note that none of the correlators shown in 
Figs. ||,|] exhibits a "volcanic" (exponential or power law 
rather than Gaussian) shape as seen for the longitudi- 
nal correlators of Refs. [0. We have not been able to 
trace the origin of this discrepancy; it may be due to the 
different source (hydrodynamics with freeze-out along a 
sharp hypersurface) used in Refs. @], but why this should 
manifest itself in this way is not obvious. From general 
arguments we would expect at small q a Gaussian be- 
haviour with a curvature related to the longitudinal size 
of the effective pion source from lu decays; the longitu- 
dinal correlators in Refs. || seem to decay much more 
steeply for small q. We have checked our results with 
two independent programs, based on the formulae given 
in the Appendices. 



VI. EXTRACTING HBT-RADII FROM THE 
CORRELATOR 

Looking at Figs. || and [| it is clear that more quanti- 
tative methods are needed to characterize the shape of 
the correlator. For an interpretation of the correlator 
in terms of the space-time structure of the source rela- 
tively small changes of its shape and its pair momentum 
dependence play an important role. One would there- 
fore like to describe the key features of C(q, K) by a 
small number of fit parameters which are sensitive to 
this space-time structure. The usual procedure is to per- 

As 



form a Gaussian fit with the functions (2.6) or (2 



we will see this method runs into systematic problems if 
the 2-particle correlator does not have a perfect Gaussian 
shape, e.g. due to long-lived resonances. Not only do the 
functions (2.6) or (2.8) fail to give a good fit, but by not 



correctly accounting for the non-Gaussian features one 
throws away important space-time information contained 
in the resonance decay contributions to the correlator. 

In this Section we discuss several different Gaussian 
fitting procedures which clearly demonstrate these dif- 
ficulties. The main reason for presenting this basically 
flawed approach is (i) that it is the method mainly used 
so far in the experimental analysis and (ii) that the dis- 
cussion throws some light on how one should compare 
HBT radius parameters extracted by different groups us- 
ing different procedures. After having understood the 
problems and the systematic uncertainties they generate 
we will then suggest a more reliable approach in the next 
Section which also accounts for non-Gaussian features in 
a quantitative way. 
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A. Two-dimensional Gaussian fits to the correlator 

We start by discussing 2-dimensional fits to C(q, K) 
with two parameters \i{K), Ri{K) [i — o,s,l). We ap- 
proximate the numerical function in the directions qi as 
follows: 

C( qi , q k ^i = 0;K) « 1 + A 4 (K) e - R ^^ , i = o,s,l. 

(6.1) 

The optimal parameters Aj(K) and Ri(K) are deter- 
mined by minimizing the following expression: 

n 

J2 C(q»,q k9ti = 0; K) - In A, + R 2 M?f = mm . 



Unidirectional Gaussian fits, Y = 0, n f = 
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FIG. 5. The Cartesian HBT radii Ri, i — o,s,l and 

their corresponding intercept parameters A;, extracted from 

the correlator C(q, K) via 2-dimensional fits according to 

Eq. (^^). Shown are results at Y = as function of K±_ for 

■qf = 0. Top row: sideward direction. Middle row: outward 

direction. Bottom row: longitudinal direction. The different 

lines indicate the effects of including various sets of resonances 

as described in Fig. [| 

The label v runs over a set of n equidistant values q\ 
between and 50 MeV for which the correlator C {q" , K) 
was calculated numerically. Although the procedure ( |6.2| ) 
is conceptually different from an experimental fitting pro- 
cedure in that the function to be fitted is known exactly 
and the resulting optimal fit parameters thus don't have 



statistical error bars, they can still vary systematically 
depending on the selection of the fit points q\ and the 
minimization function (6.2). These systematic variations 
reflect the possible non-Gaussian features of the correla- 
tor, but not in a way that allows to easily quantify them. 
As long as the deviations from Gaussian bahavior are 
small, the extracted Gaussian fit parameters i?i(K) and 
Ai(K) are expected to be useful for a simple characteri- 
zation of the main features of the correlator. 



1. No transverse flow 

For the case r\t = the results from independent 2- 
dimensional fits to the correlator in the "side" (top), 
"out" (middle), and "long" (bottom) directions are 
shown in Fig. ||. The left column shows the Cartesian 
HBT radii, the right column the associated intercept pa- 
rameters resulting from the fit, both as functions of K± 
at Y = 0. 

The fitted intercept parameters follow roughly the be- 
haviour expected from Eq. ( [3.12] ) and Fig. [l| Upon 
closer inspection one sees, however, that also some of 
the shorter-lived resonances, in particular the inclusion 
of the u>, have a significant lowering effect on A. These 
effects are different in the three Cartesian directions and 
strongest in the longitudinal direction, where even with- 
out any resonance effects A < 1 at small K±. 

The deviations of the intercept parameter from unity 
reflect non- Gaussian features of the correlator. For short- 
lived resonances these are weak, except in the longitudi- 
nal (^-direction where the correlator has been known to 
show at small K± a somewhat steeper than Gaussian fall- 
off due to the rapid boost-invariant longitudinal expan- 
sion of the source fll| l , even in the absence of resonance 
decays. The main non-Gaussian effects come from the w 
and, of course, from the long-lived resonance. The latter 
affect, however, only A and not the HBT radii extracted 
from the Gaussian fit, while the lu also changes the radius 
parameters. 

The fit accomodates these non-Gaussian fea tures by 
lowering the intercept A. 



As discussed in Sec. |VC| the 
main origin of non-Gaussian effects due to resonances is 
the tail in the time-distribution of the decay pions. Ac- 
cording to Eqs. (2.7) this is expected to affect R a and Ri, 
but not R s . Eq. (5.3) tells us, however, that the "out" 
and "side" -behaviour of the parent resonance distribu- 
tion gets mixed in the pair distribution of the daughter 
pions, so some fraction of this effect propagates into the 
side-correlator of the decay pions. On the other hand 
there remains the fact that, compared to R s , in R a an ad- 
ditional lifetime effect comes in through the term fi\(t 2 ) 
in (2.7b); this contribution increases quadratically for 
small values of K±, saturating above K± = m n where 
Px ~ 1- This explains very nicely the initial drop and 
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subsequent rise of A in the outward direction, which is 
particularly prominent for the u> contribution. 

Let us now turn our attention to the HBT radii in the 
left column of Fig. |5| and begin with a discussion of R s . 
Its size remains essentially unaffected by the short-lived 
resonances with lifetimes of order 1 fm/c, but the u> af- 
fects R s . This effect dies out rapidly for increasing M± 
due to the decreasing w-fraction f u (K± , 0) , but the re- 
sulting M±- dependence of R s complicates the extraction 
of the transverse flow from it jyj|l3| • 

The origin of the effect has already been qualitatively 



explained in Sec. VC and above by referring to Eq. (5.3) 



A somewhat more quantitative estimate can be obt ained 
by studying the space-time variances of Sec. 



Ill 



even 

though the discussion presented there makes it clear 
that this will provide only an upper estimate for the u>- 
contribution to R s . Considering only the direct pions and 
those from w-decays and calculating R 2 S = (y 2 ) according 
to (3^8) we find 



(y 2 > = /dir(y 2 



Uv 2 



(6.3) 



with 



<2/ 2 L = 



/ d'x E ± Jh Jo 00 drT e^V S^lx - P 



E± T p± 



f R fd*xS^(x,P±) 



Using f dir + f u = 1 this yields 



R ~t~ /corr fu 



l 

fa • 



(6.4a) 
(6.4b) 

(6.5) 



This result can be explained as the effect of the oj prop- 
agating in the y-direction before decaying or, more for- 
mally, as the effect of the "out" -"side" -mixing in the de- 
cay kinematics expressed by Eq. (5.3). Numerically, we 
determined the factor / corr « 0.52 at K± = which 
leads to / corr • f u « 0.1 at the same point. Putting this 
together with the width of the w-resonance, p = 23.4 fm, 
one obtains for the side variance W (y 2 ) = 8.9 fm. 

This is obviously much larger than the 5.5 fm extracted 
at K± = 0, since the curvature (2.11) does not coin- 



cide with the fitted width. For longer living resonances 
this discrepancy will, of course, be even larger. Another 
number to compare with is the half width i^ alf of the 
correlator C(q s ) at q o =qi=J£=0, including all short-lived 
resonances plus the u>. We find the hierarchy 



8.9 fm > R 



half 



6.4 fm > R s w 5.5 fm . (6.6) 



We conclude that estimates of resonance effects, based 
on space-time variances like yj (y 2 ), as e.g. done in |L8|, 



are quantitatively unreliable. The half width i?^ alf is 
close to the result one would obtain from a Gaussian 
fit to the correlator when the intercept A is fixed to the 
value of ( |3.12 ) (as e.g. done, albeit simultaneously in all 
three g-directions, in Refs. dj4l|]). The difference to our 
procedure which lets A float is significant, and since C 
at q = is not experimentally accessible, a comparison 
of i?^ alf with data El| is clearly dangerous. The authors 
of Refs. H also find that at low K± resonance decays 
can increase the longitudinal HBT radius Ri by up to a 
factor 2; in a Gaussian fit with floating A we never see 
resonance induced increases in Ri by more than 1.5 fm. 
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FIG. 6. The correlator C(q 3 ) at q o =qi=K.=0, taking into 
account only direct pions and pions from uj decays. The upper 
panel shows the three contributions according to Eq. (3.£ ): di- 
rect-direct pairs (dashed), direct-direct + direct-u; (dotted), 
and all contributions (including the u — ui term where both 
pions come from uo decays) (solid line). The lower panel com- 
pares the same solid line to different Gaussians whose radius 
parameters correspond to the curvature at q a — (dotted 
line), the half width of C(q a ) (dash-dotted line), and the op- 
timal Gaussian fit according to Eq. ( |a.2| ) (dashed line). 

We have compared the Gaussians corresponding to the 
numbers given in Eq. (|6.6|) with the true "side" -correlator 
in the lower panel of Fig. [| In the upper panel we show 
t he th ree contributions to the correlator (see (3.3) and 
(JTlQh) coming from pairs of two directly emitted pions 
(dashed line) , from pairs of one direct and one oj decay 
pion (difference between dotted and dashed lines), and 
from pairs where both pions come from to decays (differ- 
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ence between solid and dotted lines). While it is obvious 
that the u> contributions are concentrated at lower q s - 
values than the direct one, the tail from the mixed direct- 
u) contribution is still appreciable outside the half-point 
of the direct term near q s = 30 McV. It therefore ap- 
pears impossible to cleanly separate the correlation func- 
tion into "core" and "halo" contributions with different 
g-support [i5fl . In particular, the recent suggestion by 
Csorgo ]|lj to extract the "core" radius by performing 
a Gaussian fit to the g s -tail of the correlator, excluding 
the range q s < q cu t where g cu t ~ 30 MeV, is likely to run 
into systematic problems. 

We now turn to R . At K± = the two trans- 
verse radius parameters R s and R are equal by symme- 
try and all above considerations carry over to the 
"out" -direction. At non-zero K±, R Q receives an addi- 
tional contribution from the source lifetime as indicated 



by Eq. ( 2.7b ). Although for the u the use of this ex- 
pression is no longer quantitatively reliable, it gives the 
correct tendency. Short-lived resonances do not destroy 
the Gauss ian shape o f th e out-correlator, and for them 
Eq. (2.7b) (with Eq. fl3.8|)) can be used without restric- 
tions. It is obvious that even the short-lived resonances 
contribute through their lifetime to the term @\(t 2 ) in 



(2.7b), strengthening the rise of R a in Fig. || at small 
K±. 

The strongest resonance effect is seen for Ri, which is 
affected even by the short-lived resonances. These effects 
disappear for large K± due to the decreasing resonance 
fractions f r , but at small K± they are significant. Due to 
the existence of (weak) non-Gaussian features already in 
the absence of resonances the space-time variances are of 
limited use for a quantitative discussion of the effects, and 
we leave the reader with the numerical results sho wn in 
Fig. |^. A qualitative interpretation was given in Sec. V C . 



2. Non-zero transverse flow 

In Fig. ^ we show t he p ar ame ters A.;, Ri obtained from 
the 2-dimensional fit (5.1), (6.2) for the case of non-zero 
transverse expansion with ijf = 0.3. Comparing with 
Fig. H one sees that the effect of the resonances on the 
radius on the HBT radius parameters are weaker, and 
that correspondingly the non-Gaussian effects caused by 
the short-lived resonances and the lo (which lead to devi- 
ations of the intercept parameter A from unity) are less 
pronounced. In fact, the only remaining effects of these 
resonances on the HBT radii come from the terms ~ (P) 



in R and Ri (see (2.7)) and are due to the additional 
contribution to the particle emission duration from the 
resonance lifetimes. The geometrical effect of resonance 
propagation away from the direct source, described by 
the second term in (6.5), has disappeared, even for the oj. 
The reason was already discussed in Sec. VC, Eq. (5.9): 



due to transverse flow the transverse size of the emis- 
sion region for heavy resonances is smaller than that of 
the direct pions, and since they don't live very long they 
usually don't make it outside the source of direct pions 
before decaying. Thus they don't lead to an increase of 
the spatial source size. 



Unidirectional Gaussian fits, Y = 0, r| f = 0.3 



I 4 

c/) 

2 




cf" 

2 




thermal 

+ P . , 
+ A, K, L 
+ to 



6 







FIG. 

T)f = 0. 



200 400 600 200 400 600 800 

KJMeV) KJMeV) 

7. Same as Fig. H, but for non-zero transverse 



1.0 

0.8 

0.6 

0.4 S 

0.2 



1.0 
0.8 
0.6 

o> 

0.2 



1.0 



0.6 
0.4 ' 
0.2 





flow 



As a consequence, the decrease of R s with increasing 
K±, which is characteristic for transverse collective flow 
of the source pT| , ^3[ , is no longer modified by the pions 
from to decays. This is, of course, the ideal situation 
one might hope for in order to extract quantitative dy- 
namical information from HBT data. Unfortunately, the 
problem remains that, if the measurement finds a (not 
too strong) Mj_ -dependence of R s , it could still be due 
to either weak transverse flow without resonance contam- 
inations as in Fig. [?], or to w-decay contributions in the 
absence of transverse flow as in Fig. ||. (Other mecha- 
nisms like transverse temperature gradients might also 
create an Af^-dependence We must find a more 

quantitative analysis tool which allows us to tell whether 
the M^-dependence of R s is associated with w-decays or 
not. 



B. Five-dimensional Cartesian Gaussian fits to the 
correlator 

Before approaching this task in the next Section, we 
will now also discuss some generic features of multi- 
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dimensional Gaussian fits to the exact correlator where 
all HBT parameters and the correlation strength are de- 
termined simultaneously. This is clearly desirable in or- 
der to avoid the problem of having three different cor- 
relation strengths in the three Cartesian directions, as 
happens when the three radii R s , R , and Ri are deter- 



+2R 2 



(6.7) 



mined by separate 2-dimensional fits according to (3.1), 
because such a result is clearly unphysical. It is also nec- 
essary for the determination of the cross-term R D i and 
for a fit with the YKP parametrization (2 



Multidirectional Gaussian fit, Y = 
Ti f = Ti f = 0.3 
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FIG. 8. The Cartesian HBT radii R a , R 3 
cept A, obtained from the 5-dimensional fit (6.7) to G(q, K), 
as functions of K± at Y — 0. R i is not shown since at Y — 
it vanishes due to symmetry. Left column: no transverse flow, 
r/f — 0. Right column: 77/ = 0.3. 

In this Subsection we extract the Cartesian parameters 
R Q , R s , Ri , R i and A from a 5-dimensional fit which 
minimizes the expression 



E 

i/=i 



[lnG(q^,K) - In A + R 2 a (q» ) 2 + R 2 (q, 



v\2 



The label v again runs over n fit points q y which were 
chosen to lie at equal distances between and 50 MeV 
along the three Cartesian g-axes and along the two diag- 
onals (q s = 0,q o = qi) and (q s — 0,q o = -qi). This 
is, of course, different from a typical experimental q- 
distribution. They were selected to economize in the 
number of fit points where the exact correlator had to 
be computed. 

The results of the fit (6.7) are shown in Fig. |^, again 
for rjf = and r/f = 0.3 at midrapidity Y = 0. Let us 
first look at the incercept parameter A. Comparing with 
Figs. H @we see that the A- value from the 5-dimensional 
fit lies somewhere between the three different values ob- 
tained in the 2-dimensional fits. As before it reflects 
the deviations of the correlator from a Gaussian. Since 
such deviations exist in the (^-direction even without res- 
onances, due to strong longitudinal expansion, A slightly 
deviates from 1 even in the absence of resonance decays. 

The need for the fit to compromise on a unique inter- 
cept parameter affects the optimum values for the HBT 
radius parameters. For a fixed correlation function, a 
decrease of A leads automatically also to a smaller Gaus- 
sian radius as found by the fit. Since the compromise 
value for A lies above the value A;, but below the values 
A s and A Q from the 2-d fits, Ri increases and R s , R de- 
crease in the 5-d fit relative to the 2-d fit values. (This 
effect is hardly visible if resonance decays are switched 
off but becomes stronger as the resonance contributions 
(with their non-Gaussian effects) are added.) The net 
result is that even in the absence of transverse flow now 
the resonance effects on R s and on its M^-dependence 
appear quite weak. Even the resonance contribution to 
the lifetime effect in R Q (the quadratic rise at small K±) 
becomes less pronounced. 

For completeness we show in Fig. || also results at for- 
ward pair rapidity Y = 1.5. In this case the fit gives, 
of course, a non-vanishing cross-term R a i with the ex- 
pected i^-behaviour pffi] . It is affected by resonances 
essentially at the same level as R . The only other qual- 
itative difference is the much smaller value of Ri relative 
to Y — 0; this is an effect of Lorentz contraction. 

It must be stressed that the differences between this 
Subsection and the previous one are purely due to fit 
systematics. Depending on how the exact correlator is 
fitted to a Gaussian the extracted Gaussian radii show 
significant differences. Since in the experiments the in- 
tercept parameter cannot be directly measured and must 
be fitted simultaneously with the HBT radii, we adopted 
the same procedure and let A float in the fit. Schlei [^],[ll| , 
on the other hand, in his Gaussian fits has always fixed A 
at the value given by Eq. (3.12). In the presence of non- 
Gaussian effects due to resonance decays our fits give 
smaller A's and, therefore, smaller HBT radii than his 
fit. This explains why the resonance effects on the trans- 
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verse radii and their if^-dependence were found to be 
much stronger in Refs. §,|| than in our work here. Ac- 
cording to the second inequality in ( |6.6| ) the difference in 
R s is about 1 fin, in good agreement with his compared 
to our results. 

Multidirectional Gaussian fit, Y = 1 .5 
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FIG. 9. Same as Fig. [], but for forward rapidity Y = 1.5. 
Now also R i is non-zero. 



C. Five-dimensional Gaussian YKP fits to the 
correlator 

The extraction of the Yano-Koonin velocity from a fit 



according to Eq. (2.8) is a non-linear problem. To main- 
tain the simplicity of a least-square fit with linear fit pa- 
rameters we have reformulated the YKP fit problem as 
follows: We rewrite ( ]2.Sj ) in the form 

C(K, q) = 1 + A YKP (K) exp [-R 2 ± (K) q\ - ^ 3 (K)gf 
-R 2 00 (K)(q f + 2R 2 , 3 (K)q°q l ] , (6.8) 

with 







i- v 1 
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+ V 2 
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- i? R 2 



(AD 2 )' 
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1 + v 2 

^03 



p2 

-"-00 



^33 



(6.9a) 
(6.9b) 
(6.9c) 
(6.9d) 



We then proceed as with the Cartesian parametrization 
in Sec. VI B, using the same set of fit points as before, but 
expressing them through their components q±,qi, and 

R±, R33, Roo, and i?03- 



q°, in order to determine A 



Finally we solve Eqs. (6.9) for the YKP parameters. 



Multidirectional YKP fit, Y = 
r| f = n f = 0.3 



I 4 



2 


6 

i 4 

of 

2 

6 

14 
o 
rr 

2 




. 0.6 

A Y KP 

0.4 
0.2 










thermal 

" + P , , 

■ ^ v 
+ CO 

























200 400 600 200 400 600 800 
KjJMeV) K ± (MeV) 

FIG. 10. From top to bottom: the YKP fit parameters 

R±, Ro, R\\, and A YKP as functions of K± at mid rapidity 

Y = 0. Left column: no transverse flow, rjf — 0. Right 

column: rjf = 0.3. At Y — the Yano-Koonin velocity v 

vanishes exactly. 

However, the one-to-one correspondence between the 
YKP and Cartesian radius parameters does not imply 
that in a fit to experimental data both sets of fit param- 
eters can be determined with similar accuracy. At mid 
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rapidity, for instance, where q° — (3±q , the YKP fit be- 
comes for small transverse pair momentum increasingly 
insensitive to Ro, since q° — > for (3± — > 0. As a result, 
in the space of YKP fit parameters the confidence region 
for one standard deviation is very elongated in Rq. The 
actual fit value of i? thus develops a strong sensitivity 
to relatively small systematic deviations of the correla- 
tor from a Gaussian shape. Since the procedure (6.7) 
adopted here does not allow to associate errors to the 
extracted fit values, we present in Figs. |l^ and |ll] the re- 
sults only for sufficiently large values of Kj_ where such 
systematic effects were found to be small. 

For Y = (/?; = 0, Fig. [l(]) the systematic uncertainty 
at small values of K± (J3±) affects only Rq, according 
to the arguments presented above. For K± < 100 MeV, 
we found that the Rq value extracted from the Gaussian 
fit develops a strong dependence on the choice of the fit 
points while this problem disappears at larger values 
of K±. For Y = 1.5, (/?; ^ 0, Fig. [ll]), similar systematic 
uncertainties at small K± affect also i?y and Yyk- Ac- 
cordingly, the corresponding curves in Fig. [n] have been 
cut off at small K± . 

The intercept parameters A YKP extracted from the 
fit to (^^) essentially coincide with those from the 5- 
dimensional Cartesian fit. This is expected since in both 
fits the same set of fit points was used. Also, the results 
for R± compare very well with R s in the Cartesian fit. 
For a Gaussian correlator the formalism of space-time 
variances says R\ = R 2 S = (y 2 ). The equality R± = R s 
remains essentially unaffected by the non-Gaussian fea- 
tures of the correlator in the presence of resonance de- 
cays. 

The longitudinal YKP parameter Rn is affected by res- 
onance decays roughly in the same way as Ri in the 
Cartesian fit at Y — 0. This is expected because at 

Y — the two parameters are agai n ide nti cal on the level 
of space-time variances, see Eqs. ( 2.7c ), ( 2.10b ). There 
is no drastic change for Rn as one goes from Y = to 

Y = 1.5: All values (with and without resonances) de- 
crease somewhat, because one approaches the forward 
end of the source, and the longitudinal homogeneity re- 
gion thus shrinks a bit. 

The most significant resonance contribution is seen in 
the lifetime parameter Rq. This agrees with our argu- 
ments that the dominant effect from resonances on the 
correlation function arises from their finite lifetime. 

At Y = the Yano-Koonin velocity v vanishes Jl2],[l3| . 
This is reproduced by the fit. At forward rapidity v is 
non-zero. In the fourth row of Fig. ^ we plot the Yano- 
Koonin rapidity Y YK = | In fl^J as a function of the 
transverse pair momentum. For longitudinally boost in- 
variant sources, the YK rapidity is known to coincide 
with the pair rapidity, Y YK (K±,Y) = Y. For the class 



butions gave a linear relation between the two quanti- 
ties, Y^ K (K X ,Y) = c(K±)Y. The independence Y YK 
provides direct experimental access to the longitudinal 
expansion of the source. For thermalized models the 
proportionality constant c(K±) slowly approaches unity 
from below as K± increases 13 . This is clearly seen in 
Fig. [ll] which also shows that resonance decay contribu- 
tions have a negligible influence on this relation. 

Multidirectional YKP fit, Y=1.5 
T] f = T] f = 0.3 
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FIG. 11. Same as Fig. |l^, but for forward rapidity Y = 1.5. 
The additional fourth row now shows the Yano-Koonin rapid- 
ity at function of K± . 



VII. q- VARIANCES OF THE CORRELATOR 

We have seen that resonances, in particular the w 
with its intermediate lifetime, create appreciable non- 
Gaussian effects in the two-pion correlator, and that 
these deviations from a Gaussian shape can thus contain 
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additional information about the space-time distribution 
of the source and its physical origin. They also have neg- 
ative effects on the extraction of HBT radius parameters 
from Gaussian fits and affect their i^-dependence in a 
way which, within the framework of Gaussian fits, is dif- 
ficult to quantify and to control systematically. 

In this Section we therefore study an alternative ap- 
proach. We suggest to extract the HBT radius parame- 
ters and quantify the deviations from Gaussian behaviour 
by studying the normalized second and fourth order q- 
moments of the correlator C(q, K). We first develop the 
necessary formalism and then apply it to the correlation 
functions calculated from our class of source models. 



A. General formalism 

According to Sec. |l[ the most general Gaussian ansatz 
for the correlator is 



C(q, K) = 1 + A(K) exp 



i,3 = l 



D ij {K)q j 



> (7-1) 



where the qi are the three independent relative momen- 
tum components obtained after resolving the on-shell 
constraint q° = f3 ■ q. For such a Gaussian correlator, 
the HBT parameters Dy(K) can be obtained by either 
fitting the various widths of the correlator as done in 
Sec. VI, or by computing the integrals 



((QiQj)) 



Jd 3 qq iqj [C(q,K)-l} _ 1 



/d3g[C(q,K)-l] 



- (zr^K)). 



(7.2) 



and inverting the resulting matrix of second order q- 
moments. 

For a non-Gaussian correlator we may define the HBT 
radius parameters in terms of these "q- variances" : hav- 
ing determined the matrix -D(K) by inverting the matrix 
((q (g> g))(K) of q- variances, we define 
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when q s , g G , qi are used as independent coordinates, and 
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(7.4) 



if one uses instead q±,qi,q° as independent variables. 
Eq. (7.3) corresponds to the Cartesian parametrization 
(2.6), generalized to systems without azimuthal symme- 
try by allowing for non-vanishing "side-out" and "side- 
long" cross terms. Eq. (7.4) corresponds to the YKP 



parametrization (2.8) which applies only to azimuthally 
symmetric systems, and the zeroes in the matrices on the 
left and right hand side reflect this symmetry. 

Similarly, the intercept parameter can be defined in 
terms of the g-variances and the zeroth order g-moment 
as 



A(K) = 7T- 3 / 2 v'detD(K) J d 3 q [C(q, K) - 1] 



(7.5) 



which reproduces the correct value for Gaussian correla- 
tors of type (frTl]). 

The deviations from Gaussian behaviour in the correla- 
tor are then related to higher order g-moments. A general 
discussion, including their derivation from a generating 
functional from which the full correlator can be recon- 
structed, is given in Ref . |Q . Since C (q, K) is symmet- 
ric with respect to interchange of the particle momenta 
Pi and P2 and therefore even under q — > — q, all odd 
g-moments vanish. The first non-Gaussian contributions 
thus show up in the fourth order moments. 

Application of the method of g-moments thus generally 
requires at least an inversion of the matrix ( |7.2| ) for the 
determination of the HBT radius parameters and a dis- 
cussion of the 4-dimensional tensor of fourth order mo- 
ments for the non-Gaussian aspects. A complete such 
analysis in three-dimensional g-space will be postponed 
to a future publication. Here we will perform a unidirec- 
tional analysis, where these technical complications do 
not arise, and compare it with the unidirectional Gaus- 



sian fits of Sec. VI A 



We thus consider the correlator along one of the three 
axes qi (i — s,o,l or i =_L, I, 0) which we denote by C(qi), 
suppressing for simplicity the K-dcpcndcncc: 



C(q i ) = C(q i ,q j ^ = 0,K) 



(7.6) 



The HBT radius parameter in direction i and the corre- 
sponding intercept parameter are then defined as 



R 



1 



2«a 



a ,, _ fdq iq ?lC(qi)-l] 



Jd qi [C( qi )-l} 
Xi = (Ri/V^) I dqi [C( qi ) - 1] . 



(7.7a) 
(7.7b) 
(7.7c) 



To extract the moments ((gf )) from data one replaces 
( 7.7b| ) by a ratio of sums over bins in the gi-direction. 
The higher the order n of the q- moment, the more sensi- 
tive are the extracted values to statistical and systematic 
uncertainties in the region of large qi. First investiga- 
tions with event samples generated by the VENUS event 
generator indicate that the current precision of the data 
in the Pb-beam experiments at the CERN SPS permits 
to determine the second and fourth order g-moments. 
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Accordingly, we restrict our discussion of non-Gaussian 
features to the "kurtosis" 



A; 



«<z, 4 » 



1 



(7- 



3((<Z, 2 )) 2 

In the following Sec. [VII B we will study the in- 
dependence of the HBT radius parameters, the intercept 
and the kurtosis as defined by Eqs. (7.7) and (VM). 



B. Unidirectional results for the q- moments 

In this Subsection we present a numerical analysis of 
the correlation functions computed in Sec. ^ in terms 
of their g-moments along the three Cartesian directions, 
and give a comparison with the unidirectional Gaussian 



fits presented in Sec. VIA . 

|l2] shows the HBT radii ( 7.7a ) and the kurtosis 



Fig. 



( |7.8[ ) along the "side"-, "out'-, and "long"-axes (from top 
to bottom). The left and right column of plots corre- 
spond to zero and non-zero transverse flow of the source, 
respectively. In each panel we plot as the upper set of 
curves the HBT radius parameter Ri in fm, with differ- 
ent line symbols denoting the effects of including various 
sets of resonances as before. They should be compared 
with the lines shown in the left columns of Figs. ^ and 
fj], respectively. The lower set of lines (clustered around 
values near 0) denote the corresponding kurtosis Aj in di- 
mensionless units. These contain the lowest order infor- 
mation on the non-Gaussian features of the numerically 
computed correlation function. 

The comparison of the HBT radius parameters defined 
via the g-variances of the correlator with those from the 
Gaussian fit ( |S.2| ) shows a remarkable agreement. As 
stressed above, in the presence of non-Gaussian features 
in the correlator, the only well-defined defini tion of the 
HBT radii is provided by the q- variances ( 7.7a ), while the 
Gaussian fit results have possibly severe systematic un- 
certainties related to the details of the fit procedure. The 
agreement between the corresponding curves in Figs. |^, 
fj], and |l2| indicates that w e wer e "lucky" with our choice 
of fit prescription in Sec. VIA. An essential reason for 
the good agreement was our decision to let the inter- 
cept parameter A float in Eq. (6.2), i.e. to perform a 2- 
dimensional rather than a 1-dimensional fit as in Refs. Q . 
The discrepancy between the HBT radii shown in those 
papers and those shown in Fig. ^ thus simply reflect the 
systematic uncertainties of extracting a Gaussian width 
parameter from a non-Gaussian correlator. In view of 
these uncertainties, the existence of a clear-cut definition 
via the q- variance of the correlator becomes crucial. 

The space-time interpretation of the HBT radius pa- 
rameters has so far been largely based on their relations 



(2.7), (2.10) with the space-time variances of the source 
which are only true for Gaussian correlators. The agree- 
ment between the HBT radii from g-variances and from 



(appropriate) Gaussian fits suggests that these relations 
continue to be useful for the space-time interpretation of 
the correlation functions. 

Unidirectional q-variances, Y = 
il f = n f = 0.3 
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FIG. 12. 1-dimensional q- variances according to Sec. VII A 
The plots show for the three Cartesian directio ns i = s,o,l 
the radius parameters Ri(K±) defined by Eq. (7.7a) (upper 
set of c urves in each panel) and the kurtosis Ai(K±) defined 



by Eq. (7.7a) (lower set of curves in each panel). The radii 
are given in fm, the kurtosis in dimensionless units on the 
same scale. Left column: no transverse flow, 77/ = 0. Right 
column: 77/ = 0.3. The pion pairs have rapidity Y — in the 
CMS. 

In view of the above agreement between the two types 
of HBT radius parameters, and of our discussion of the 
interplay between the values of Xi and Ri in various types 
of Gaussian fits to a given correlation function, it is not 
surprising that the intercept parameters extracted from 
Eq. (fr?7cl) also agree very well with the ones extracted 
from the unidirectional Gaussian fits and shown in the 
right columns of Figs. || and |?]. They are therefore not 
presented again. 

The interesting new information is, of course, con- 
tained in the kurtosis Aj and their .KJ-dependence shown 
in Fig. [l^. In the side direction the appearance of a non- 
vanishing (positive) kurtosis is clearly linked to the influ- 
ence of the w decays on the correlation function and to its 
visibility in the HBT radius parameter R s . This implies 
that the question whether or not a given isT^-dependence 
of R s is caused by resonance decays or not can be eas- 
ily answered by checking the kurtosis of the correlation 
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function. If the kurtosis vanishes (or is slightly negative) , 
it is not the to which causes the K j_-dependence. At least 
for the model studied here, the kurtosis provides thus the 
cleanest distinction between scenarios with and without 
transverse flow. Its value and if ^-dependence are thus 
very important ingredients for the interpretation of 2- 
particle correlations. 

The situation is slightly more complicated in the out- 
ward direction: as long as the source does not expand 
transversally (77/ = 0), the visibility of resonance decay 
effects in R a is clearly linked to a non-zero positive kurto- 
sis of the correlator, and vice versa. For non-zero trans- 
verse flow, however, the outward correlator begins to de- 
velop small deviations from a Gaussian [fill even without 
resonance decays; these show up in a negative value for 
the kurtosis. This effect increases for larger transverse 
pair momenta Kj_. 

The kurtosis generated by collective expansion is par- 
ticularly prominent in the longitudinal direction where 
flow-induced non-Gaussian features they have been no- 
ticed first jy]]. The bottom row of Fig. [l2| clearly shows 
the interplay of non-Gaussian features induced by res- 
onance decays (leading to a positive kurtosis) and lon- 
gitudinal expansion flow (causing a negative kurtosis). 
At small K± the resonance contributions dominate; at 
large K± the resonances loose importance while the flow- 
induced kurtosis becomes stronger, leading to overall neg- 
ative values of the kurtosis. 



VIII. CONCLUSIONS 

Within a broad class of model emission functions for 
locally thermalized and collectively expanding sources we 
have presented a comprehensive study of resonance de- 
cay effects on two-pion Bose-Einstein correlations. We 
have found that, with regard to their influence on the 
correlation function, the resonances can be subdivided 
into three classes: 

• Long-lived resonances with width < 1 MeV can not 
be resolved in the correlation measurement; they 
reduce the correlation strength A but otherwise do 
not influence the shape of the correlation function 
in the region where it can be measured. 

• Short-lived resonances with width > 30 MeV: they 
decay into pions close to their production point and 
thus do not change the spatial width of the pion 
emission function. Hence they do not affect the 
sideward correlator whose width is defined by the 
transverse spatial size of the source. In the outward 
and longitudinal correlator and in the lifetime pa- 
rameter Rq of the YKP paramctrization, which arc 
all in one way or other sensitive to the lifetime of 
the source, they contribute via the additional time 



duration of pion emission due to their own lifetime. 
These contributions are small and on the order of 
the resonance lifetime. 

• The to meson. With its width of about 8 MeV 
it is not sufficiently long-lived to escape detection 
in the correlator, but also not sufficiently short- 
lived to not change the spatial width of the emission 
function. As a consequence it can lead to severe 
non-Gaussian distortions of the correlator. 

These latter distortions cause serious problems. We 
have shown in Sec. [Vl] that both the method of extract- 
ing width parameters from the correlator via Gaussian 
fits and the calculation of these parameters in terms of 
space-time variances can lead to quantitatively unreliable 
results. The systematic uncertainties of Gaussian fits to 
non-Gaussian correlators were identified in Sec. VI (c.f. 



the discussion following Eq. (6J3)) as the primary reason 
for previous claims of much larger resonance effects on 
the two-pion HBT radii than found by us. To remove the 
ambiguities associated with non-Gaussian features of the 
correlator we have introduced in Sec. VII an alternative 



definition of the HBT size parameters and of the inter- 
cept parameter A in terms of q-moments of the correlator 
which does not rely on the assumption of a Gaussian cor- 
relator. For sufficiently high statistics data, HBT radius 
parameters determined in this way are free of systematic 
uncertainties. For the examples studied here, they show 
a much weaker influence from resonance decays than we 
had expected on the basis of previous work @]. 

The normalized fourth order g-cumulant (kurtosis) 
serves as a quantitative lowest-order measure for the non- 
Gaussian features of the correlator. It is sensitive to both 
resonance decays and flow which (at least for the models 
studied here) contribute, however, with different signs. 
The kurtosis thus provides the cleanest signal to dis- 
tinguish between scenarios with and without transverse 
flow. 

Our detailed numerical model study of g-moments has 
shown that resonance decays which modify the HBT ra- 
dius parameters (defined via the g-variance of the corre- 
lator) also lead to a positive kurtosis. It can be related 
to the long non-Gaussian tails in the source distribution 
generated by the decay pions. Collective expansion, on 
the other hand, generates a negative kurtosis because it 
tends (in our model) to let the source at its edges de- 
cay more steeply than a Gaussian. We see practically no 
flow effects on the kurtosis in the sideward direction, a 
weak effect due to transverse expansion in the outward 
direction, and a somewhat larger effect due to the strong 
longitudinal expansion in the longitudinal direction. In 
the transverse direction resonance effects on the HBT ra- 
dius R s can thus be directly correlated with a non-zero, 
positive kurtosis. The existence or not of a non- vanishing 
kurtosis A s and its ifj_-dependence can thus be used to 
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assess the amount of contamination in R s from w-decays 
and to separate these effects from transverse flow. 

g-moments thus provide significantly improved infor- 
mation on the shape of the correlation function in terms 
of a still small number of relevant parameters Xi, Ri, Aj, 
whose size and momentum dependence lends itself to an 
interpretation in terms of the geometric and dynamic 
space-time structure of the emitting source. They are 
thus expected to further adapt the HBT method to the 
increased demand for accuracy in view of the complicated 
nature of the dynamical sources created in relativistic 
heavy ion collisions and of the drastically improved qual- 
ity of recent correlation measurements. The new method 
has been demonstrated to work very well in theory. In 
view of the new high precision data from the Pb-beam 
at the CERN SPS, it appears to be experimentally fea- 
sible. It will be interesting to see how far the additional, 
higher order HBT observables improve our picture of the 
spatio-temporal evolution of heavy ion collisions. 
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APPENDIX A: THE EMISSION FUNCTION FOR 
RESONANCE DECAY PIONS 

Here, we give details of how to compute the emission 
function S r —,- K (x,p) for resonance decay pions from a de- 
cay channel r. We follow the treatment in |l(],[l7j with 
some notational improvements. The resonance r is emit- 
ted with momentum P at space-time point X^ and de- 
cays after a proper time r at x M = X^ + -p-r into a pion 
of momentum p and (n — 1) other decay products: 



The decay rate at proper time r is Te~ TT where T is the 
total decay width of r. Assuming unpolarized resonances 
with isotropic decay in their rest frame, S r ^.„{x,p) is 
given in terms of the direct emission function S^ 1T (X, P) 
for the resonance r by 

f s + f d 3 P 

S r ^(x;p) = M J dsg(s) j — 5(p-P-E*M) 



d 4 X / drT 



-IV 



S W 



X 



p 

M 



S™(X,P). (A2) 



Variables with a star denote their values in the resonance 
rest frame, all other variables are given in the fixed mea- 
surement frame. Here s — (X)"=2-Pi) 2 m the squared 
invariant mass of the (n — 1) unobserved decay prod- 



ucts in ( Al). It can vary between s_ = (J2i=2 mi ) ano - 
s+ = (M — m) 2 . g(s) is the decay phase space for the 
(ri — 1) unobserved particles. E* is the energy of the ob- 
served decay pion in the resonance rest frame and is a 
function of s only: 



E* = Jm 2 +p* 



P 



— m) — s 

2M 



(A3) 



We choose for the observer frame a Cartesian coordi- 
nate system in which the transverse momentum pj_ of 
the decay pion has only an x ("out") and no y ("side") 
component: 

j/ 1 = (E,p x ,p yi p L ) = (mj_coshy, p x , 0, m^sinhy) . 

(A4) 

In this coordinate system the resonance 4-momentum P 
is parametrized by 

P^(E p ,P X ,Py,P L ) 

= (Mi cosh Y, Pj_ cos $, Pi sin $, Mi sinh Y) . (A5) 



The first S- function in (A2) implements the energy- 
momentum constraint p ■ P = E*M. For p± ^ it can 
be used to fix the azimuthal angle <& of the resonance 
momentum P to 



<I> ± = ±<I> with 

EE P -p L P L 



cos<I> 



E*M 



P±P± 

m±M± cosh(F - y) - E*M 



P±P± 



(A6) 



We denote by P ± the two values of P obtained by in- 



serting the two solutions (|Aq) into (A5). Rewriting the 



7T + C 2 + C 3 



^-function as S(p-P-E*M) = £± p X" s f n ± 2 ± and doing 
( A1 ) the ^-integration in d 3 P/E p = M ± dMj_dY d$ we find 
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dY 



Y- 



M 



dMi 



M 2 



dTTe- rr S^ x-[X- 
x^ dir (X,F ± )<i> r ^(P ± ;p) 



d 4 X 
P± 

H 7 



(A7) 



where 



ds- 



Mg(s) 



\JP\v\ - [E*M - m±M ± cosh(Y - y)f 



(A8) 



is the decay probability for a resonance r with momentum 
P into a pion with momentum p. It is normalized to the 
branching ratio br^^ for the channel (Al) according to 

J dydp 2 ± <i> r ^ 7I {P;y,p 1 _) = b r ^ 7r . (A9) 

The case p± = is a little special: then the con- 
straint p-P = E*M in (A2) cannot be used to do the 
^-integration, but the Mi-integral can be done: 



S r -, w (x;y,p± = 0) 



= M 



dsg(s) 



2 71 



dY- 



ME* 



m 2 cosh (Y — y) 



d 4 X / drTe- rT 5^ 



x - [X + —T 
1 M 



xS?*(X,P) 



(A10) 



M, =- 



In the following we discuss only the case p± ^ 0. The 
kinematic limits for the integrals in (A7) and (A10) are, 
for given y, m± of the decay pion, determined by the 
zeroes of the square root in (A8): 

M 1± =I i± AMx, E*M m± cosHY-y) 



to^_ cosh (Y — y) — p 2 



±- 



Mpj_ JE* Z +p 2 x -mj_ cosh" (Y-y) 



cosh 2 (Y — y) — p\ 



Y ± =y±AY = y±\n 



P_ 



(All) 
(A12) 



With these ingredients Eq.rtA2) can be rewritten as 



Sr^ix.p) = M 



dsg(s) 



dTTe~ TT 



M 



dY 



r,P d 



i\/f 



J-2 



where the sum is over the two allowed value s (|A6| ) for $. 
Rewriting the square root with the help of ( A10| ) as 

1 1 



cosher ~y)~p\ y/(AM ± )* - (M ± - M ± y 

(A14) 

and introducing new integration variables v G [—1,1], 
C £ [— 7r, 7r] via 

Mi = Hi + AMi cos C , (A15) 
Y = y + vAY, (A16) 

Eq. ( A13[ ) can be further transformed into 

S r ^(x,p) = Y, J r J o °° dTTe- rT S?" (x -^ T ,P^j , 

(A17) 

with the following shorthand for the integration over the 
resonance momenta: 

el 



= M 



R 



AYdv 



dsg{s) 

J m 2 ^ cosh 2 {v AY) — p\ 
d( (M± + AM± cos C) • (A18) 



For the calculation of the correlation function we need the 
Fouri er tra nsform of the emission function. It is obtained 
from ( A17 ) as 

S r ^(q,p) = ( d A xe^ x S r ^{x,p) 



E 



d(Tr) exp 



r Jo 



-Ft 



■ <1 



•P± 



Mr 



x / d A xe iq ' x Sf^x^^ 



E 



p ± S r (q, P ± ) . 



R 1 — i 



(A19) 



Air 



where in the first step we shifted the ^.-integration 
variable and in the second step we performed the r- 
integration. For two- and three-body decays, this reads 

• For two-body decays: 
b 



9(s) 
S r ^(q,p) 



Airp* 
Mb 
Airp* 



E 



(A20) 



AYdv 



cosh 2 (uAF) — p 2 



Qq 



AIT 



1 - iQq 

(q° coshY — qi sinhy) 



^p]_{Ml - M 2 ) - [E*M - m±M± cosh(T - y)\ 



, (A13) 



--^(? cos$ ± + g s sin$ ± ) . (A21) 
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• For three-body decays (s_ = (m.2 + m 3 ) 2 , s + 
(M-m) 2 ): 



p2n 

—i / d?/> exp ( sinh 774 cos V') 
•/ 

x sin (rgj_ cos <^ cos ip — rq± sin sin ip) 
= 7T (/ (VC-i£>) - I (VC + »!>)) , (B3b) 



Mb y/[s - (m 2 + m 3 ) 2 ][s - (to 2 - m 3 ) 2 
2tts 



Q(M, m, m 2 ,m 3 ) 



(A22) 



Q(M, m, m 2 , m 3 ) 



y/s + - s'y/s^ ~ S 'V ( m 2 - W3) 2 - S' ■ 

bM 2 



2itQ(M, to, to 2 , to 3 ) 



x / \/[s - (to 2 + m 3 ) 2 ][s - (m 2 - m 3 ) 2 ] 



AF du 



1 yro^ cosh 2 OAF) - p]_ 
Ja 1 - iQq 



APPENDIX B: THE FOURIER TRANSFORM OF 
THE EMISSION FUNCTION 

Here, we give details of the calculation of the Fourier 
transform Sf\q,P) = J d i x e iqx Sf r (x, P) for the reso- 
nance emission functions ( [4.6| ) . The r-integration can be 
done analytically: Using q ■ x = t A ~ q Q x — q s y with A 



from (5.2c) we obtain 



• cLt e lAT exp 



(r-ro) 5 



2(At) 2 

V2^(Ar) 2 e lAr ° e^^rf / + ii4 ( Ar ) 2 ) . ( B l) 



The angular integral is also easily done: writing q Q — 
q± C0S(y9, q s = q± simp, such that q x + q s y — rq± cos(0 — 
ip) (where <p is the polar angle of x and y) , the integral 
over the angle-dependent part of the source function (4.6 ) 
is written as 



2tt 



d<p e~ irq± ^(^f) sinh *7* cos(0-$) 



dtp e~ iTq± cos (V'+0) sinh)) t cosi/ 1 



(B2) 



with -0 = = vs. Separating real and imaginary 

parts one obtains modified Bessel functions f45| : 

/ d^ exp ( £j- sinh 774 cos tp) 
j 

x cos (rq± cos cos t/> — rq± sin y> sin ip) 
= 7T (/ (VC-iD) + I (VC + *£>)) , (B3a) 



where C and -D are given in ( p.2| a,b). The remaining 
integrals over r and 77 are given in ( |5.l| ) and must be 
done numerically. 

The single particle spectrum is obtained by evaluat- 
ing S(q, P) at q = 0. Then also A, A q and D van- 
ish (i.e. the dependence on the polar angle $ of the 
transverse momentum Pj_ drops out), and C reduces 
to C — (P±suihr]t(r)/T) 2 . The transverse momentum 
spectrum is obtained by additionally integrating over the 
rapidity Y associated with P. This integral can again be 
done analytically: 



dNS 



dM 2 



dYS?*(0;M ± ,Y) 

2J r + 1 f£r 

■ M±tq e T 



■ dr i 



(2tt) 3 / 2 

xl sinh 77 t (r)) J d V exp f- ^JL- ^ 

x y dF cosh(?7 - F) 

x exp [— cosh 774 cosh(?7 — y)l 



2J r + 1 
2tt 



(2r A77)eT" Mi 



• dr i 



xtf x (^ cosher)) /o sinhr? 4 (r)) . (B4) 



The i-Ti-function results from the last integral in (B4) 
after a simple shift of the integration variable, and the 
remaining Gaussian integral over 77 is trivial. 
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